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Abstract. We study a chain of infinitely many particles coupled by nonlinear 
springs, obeying the equations of motion 

(?n = V'(g n +i - q n ) - V'(q n - q n -i) 
with generic nearest-neighbour potential V. We show that this chain carries exact 
spatially periodic travelling waves whose profile is asymptotic, in a small-amlitude 
long-wave regime, to the KdV cnoidal waves. The discrete waves have three 
interesting features: (1) being exact travelling waves they keep their shape for 
infinite time, rather than just up to a timescale of order wavelength -3 suggested 
by formal asymptotic analysis, (2) unlike solitary waves they carry a nonzero 
amount of energy per particle, (3) analogous behaviour of their KdV continuum 
counterparts suggests long-time stability properties under nonlinear interaction 
with each other. Connections with the Fermi-Pasta-Ulam recurrence phenomena 
are indicated. Proofs involve an adaptation of the renormalization approach of 
1 1 2| to a periodic setting and the spectral theory of the periodic Schrodingcr 
operator with KdV cnoidal wave potential. 



Keywords: Fermi-Pasta-Ulam problem, Korteweg-de Vries equation, cnoidal wave 
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1. Introduction 

The Fermi-Pasta-Ulam model consists of a one-dimensional chain of particles coupled 
by nonlinear springs, obeying the equations of motion 

'in = V'(q n+ i - q n ) - V (q n - q n -i)- (1) 

Here q n (t) is the displacement of the n th particle out of equilibrium at time t, and V is 
an anharmonic potential such as V{<j>) = 4> 2 /2+a<p 3 /3 with a^O (the FPU-a chain) or 
V(4>) = 2 /2+/3</> 4 /4 with p 7^ (the FPU-/3 chain). This model provides a fascinating 
paradigm of nonlinear Hamiltonian many-particle dynamics. On the one hand, it 
is simple enough to allow insight by rigorous mathematical analysis. On the other 
hand, it already exhibits a rich spectrum of phenomena of wider importance: coherent 
signal- and energy transport (as seen in biomolecules such as DNA); near-integrable 
behaviour (as documented by the fact that in certains regimes, the FPU model is well 
approximated by the Korteweg-de Vries equation, of which more below); dispersive 
shocks (as seen in molecularly resolved gas dynamics); and statistical irreversibility 
and thermalization effects despite microscopic reversibility (as described by statistical 
mechanics) . 
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A central numerical phenomenon in the system ([I]) is a crossover from energy 
trapping in a few long-wave modes at low initial energy per particle to ergodic-like 
spreading of energy to short-wave modes at high initial energy per particle. See 
Fermi, Pasta, Ulam and Tsingovjf] [1TJ] for first observations at low energy, Israilev and 
Chirikov 28J for first observations at high energy, Dreyer and Herrmann [7] for energy 
transfer to short-wave modes via dispersive shocks, and [4] for a nice review. 

A significant amount of understanding of the recurrent and non-statistical 
behaviour at low energy has emerged via approximation of the FPU model by 
completely integrable systems. In this paper 

- we argue that this level of understanding, reviewed below, is not completely 
satisfactory because it fails to cover the physically most interesting regime, fixed 
nonlincarity, nonzero energy per particle, and larger and larger system size 

- we hope to convince readers that a deeper understanding could come via establishing 
existence, and long-time stability under interaction with each other, of spatially 
periodic and quasi-periodic waves in the infinite FPU chain 

- and we rigorously carry out a first step, establishing existence of certain spatially 
periodic travelling waves in FPU which are good candidates for such special interaction 
properties, namely waves which are, in a small-amplitude long-wave limit, asymptotic 
to the celebrated KdV cnoidal waves. 

In the remainder of this Introduction we review current theoretical understanding 
of non-statistical behaviour of FPU at low energy, and informally describe our results. 

1. The invariant tori explanation. One line of thought, going back to [28], is 
to consider a fixed and finite number of particles only, and linearize, i.e. approximate 
the anharmonic interaction potential by a purely harmonic one. The resulting system 
is, of course, completely integrable, with the phase space being foliated by the invariant 
tori given by the set of states with a given fixed amount of energy in each normal mode. 
By Kolmogorov-Arnol'd-Moser (KAM) theory, when switching on the anharmonic 
terms, many invariant tori survive. This prevents typical solutions from spreading 
energy to the entire phase space. This argument can be made rigorous by using 
Birkhoff normal forms |41) . 

2. The soliton explanation. Another approach, going back to Zabusky and 
Kruskal |50j . is to keep the number of particles infinite, consider a suitable small- 
amplitude long- wave regime, and observe (see Remark 2.1 below for their heuristic 
argument) that ([lj is then well-approximated by a nonlinear, infinite- dimensional 
integrable equation, the KdV equation 

V"'(0) 

u t + 12 yjyjoj uu x + u xxx = 0. (2) 

(More precisely, the latter approximates the FPU-a chain; in case of the FPU- 
j3 chain one obtains the mKdV equation.) The explanation of non-statistical 
behaviour advocated in [SU] then went as follows: the latter equation carries 'solitons', 
i.e. solitary waves which propagate exactly under the nonlinear KdV dynamics; 
numerically these solitons also show persistent shapes and velocities after nonlinearly 
interacting with each other; and general spatially localized initial data can be well 
approximated by superposition of finitely many such solitons. Hence, again, energy 
is prevented from spreading to the entire phase space. Significant aspects of this 
picture has nowadays become rigorous mathematics. In particular, the inverse 
scattering transform introduced in 1974 by Gardner et al. |17j allows to prove the 



| Mary Tsingou was involved in the numerical work as acknowledged in the original report. 
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radiationless interaction of solitons, as well as show that all spatially localized initial 
data asymptotically split up into a superposition of solitons (see e.g. [20 for more 
information and further references). And the approximation of ([I]) by ^ can be 
made rigorous up to timescales of order (wavelength) 3 for general localized solutions 
[18] . and globally in time for special solutions of soliton type O HH US] and 
multi-soliton type [38l [37] . 

From a physical point of view, neither the invariant tori explanation nor the 
soliton explanation are completely satisfactory. This is because they are limited to 
regimes satisfying certain undesirable restrictions: 

• The invariant tori explantion is based on finite-dimensional KAM theory, and 
hence the allowed size of the anharmonicity tends rapidly to zero as the particle 
number gets large |4?j . It thus does not apply to the natural situation of many 
particles interacting via a fixed nonlinear potential. Very interesting results 
are available on KAM theory for perturbations of infinite-dimensional systems 
[29], but perturbing the harmonic lattice (or indeed other infinite-dimensional 
integrable systems like the Toda lattice or KdV) into FPU lies well beyond the 
scope of these results. 

• The soliton explanation does apply to infinitely many particles interacting via a 
fixed nonlinear potential, but it is only valid for spatially localized states; but 
spatial localization forces the energy per particle to be zerojf] 

How, then, could we access the behaviour of infinite FPU chains with fixed 
nonlinear potential and nonzero energy per particle? Our own, admittedly rather 
modest, contribution towards this challenge is the following. We show that ([I]) carries 
exact spatially periodic travelling waves 

r n {t)=R{n-ct) (3) 
with small-amplitude long-wave profile 

R{n-ct) = £ 2 $(e(n - ci)) + 0(e 4 ) (4) 
and near-sonic velocity 

c 2 = U"(0)(l + ^), (5) 

where r n (t) = q n+ i(t) — q n {t) is the relative displacement between neighbouring 
particles and $ is a KdV cnoidal wave profile of speed 1, i.e. a spatially periodic 
function such that u(x 7 t) = &(x — t) solves the KdV equation ([2]). Such solutions to 
KdV exist and are known explicitly via algebraic and geometric methods, see Figure [T] 
and the end of this Introduction. 

While interesting non-explicit existence results on periodic FPU travelling waves 
have been obtained previously [46] [TTJ |39l [21] , the key point which makes the new 
result ([3])~([5]) promising in the context of energy trapping in FPU at low energy per 
particle is now the following. The discrete cnoidal waves ([3]) _ ([5]) are not arbitrary FPU 
travelling waves with finite energy per particle, but are good candidates for exhibiting 
special stability properties under nonlinearly interacting with each other, on account 
of analogous remarkable properties of their continous KdV counterparts. As made 

§ Curiously, the numerical KdV simulations in initial paper |50j on the soliton explanation were 
done at finite energy per particle, via imposing periodic boundary conditions, but the subsequent 
theoretical analysis of soliton interactions was done at zero energy per particle, via just considering 
a fixed number of localized solitons on the whole real line. 
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Figure 1. Different cnoidal wave solutions of |2), all moving with speed one. 
We have set V"( 0) = V"'(0) = 1. The three curves represent the solutions given 
explicitly by i |4"i| , with elliptic modulus k 2 = 0.3 (dotted), k 2 = 0.6 (dashed), and 
k 2 = 0.999 (solid). The corresponding values of the period are chosen to make the 
wavespeed equal to one, i.e. 2L = &K(k)(l — k 2 + fc 4 ) 1 / 4 , and are approximately 
6.5, 7.3, and 19.4. 



more precise below, the KdV cnoidal waves are the spatially periodic analoga of the 
KdV solitons, and appropriate nonlinear superpositions, the so-called finite-gap KdV 
solutions, are (spatially quasi-periodic) analoga of the KdV multi-solitons. A long- 
term goal, then, which lies well beyond the scope and tools of this article, would be to 
investigate existence and long-time stability of discrete FPU-analoga of general finite- 
gap KdV solutions. Perhaps infinite-time stability of such waves in the spirit of the 
deep recent results of Mizumachi [351 ISZ] on discrete FPU-analoga of multi-solitons 
no longer holds, but it is conceivable that long finite-time results analogous to those 
of Hoffman and Wayne [23 on FPU solitary wave interactions can be transferred to a 
quasi-periodic setting. 

The above line of thought also suggests, as to the best of our knowledge has 
not been pointed out previously, that the linear normal modes which were nonlinearly 
evolved in the original FPU simulations [TU] and underlie the "KAM explanation" [5T] 
enjoy special interaction properties under nonlinear FPU dynamics. This is because, 
in the limit of low amplitude, the KdV cnoidal waves become close to linear waves 
plus a constant; see the dotted curve in Figure [l] and (41 1 in Section [4] 

Next, let us compare our result to previous work on localized waves in the 

FPU model. Existence of travelling waves with constant asymptotic values at infinity 
has been established via different approaches: variational methods [TBI 14*51 [TT1 |3T)1 |2"T] . 
see also [32 [22] for a generalization to nonconvex potentials, center manifold 
arguments [25], and comparison to KdV [12] . The latter two approaches are limited 
to small amplitude, but also deliver the waveform. 

We close this Introduction by describing briefly the above-mentioned finite-gap 
KdV solutions, and explaining how they encode and reveal remarkable nonlinear 
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interaction properties of the KdV cnoidal waves in Figure [T] The finite-gap solutions 
can be constructed via algebraic and geometric methods going back to Its, Matveev, 
Dubrovin, McKean, and van Moerbeke [251 [S7J GO HI ISS] , and are given by 

u(M)=c + 2^1ne(0, (6) 
6(£;B)= exp{27ri(k,0 +7ri(k,Bk>}, 

where (., .) denotes the scalar product in R and c is a constant. Moreover, B is a 
symmetric matrix and £ — . . . ,£;v) denotes the vector of the generalized phases 
£j = KjX — Ldjt + 4>j, j = 1, . . . ,N, which are determined by an underlying prescribable 



Riemann surface (see Appendix B I . In general, the Kj 's and LUj 's are incommensurable 
quantities and thus KdV solutions of the form ^ are quasi-periodic in x and t. In the 
one-gap case (N = 1) the solutions Q reduce to the, spatially periodic, KdV cnoidal 
waves given by 

u{x, t) = E 2 + (E 3 -E 2 )cn 2 (J ^ ~ El (x-2(E 1 +E 2 +E 3 )t); k 2 ) , (7) 



where E\ < E 2 < E 3 and k 2 = (E 3 — E 2 )/(E 3 — Ei) is the elliptic modulus. Here 
cn denotes one of the Jacobi elliptic functions. See Figure [l] The solutions ^ were 
already known to Korteweg and de Vries in 1895 ^30 . 

Physically, the general finite-gap solutions (|6| can be interpreted as a linear 
superposition of N cnoidal waves of the form |7| plus a nonlinear interaction term. 
Here the diagonal elements of the symmetric matrix B determine the periods of these 
cnoidal waves, while the off-diagonal elements give their nonlinear interaction. In this 
precise sense, these solutions are finite-energy-per-unit-volume analoga of the KdV 
multi-solitons. 

For a more detailed description of quasi-periodic KdV solutions we refer 
to |Appendix B| 



2. Main result 



We now describe precisely or result on existence and shape of cnoidal-type waves in 
Fermi-Pasta-Ulam chains. 

The governing equations are 

d 2 t q{j,t) = V l {q{j + l,t)-q{j,t))-V l {q(j,t)-q{j-l,t)) (jeZ),(8) 

where q(j, t) denotes the displacement out of equilibrium of the j'th particle. Formally, 
the associated Hamiltonian energy 

H{t) = E ( P -^Y~ + + *) - <?(•?> t))) (9) 

(with dtq(j,t) =: p(j,t) denoting the particle momenta) is conserved along solutions, 
but our interest is in infinite-energy solutions. In this paper we make the following 
assumptions on the nearest-neighbour potential: 

v e c 4 , 7(0) = v'(o) = o, v"(o) > o, y"'(o) > o. (10) 

Well-known examples are given by the potential of the cubic FPU chain V(r) = 
ar 2 /2+br 3 /6, the Lennard-Jones-(12,6) potential V(r) = Ar~ 12 — Br~ 6 , and the Toda 
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potential V(r) — a(e~ l3r + (3r — 1). The latter gives rise to a completely integrable 
Hamiltonian system. Let us denote by 

r(j,t)=q(j + l,t)-q{j,t) 

the distortion of the jth bond length out of equilibrium. Moreover, we introduce the 
shift operators S ± q(j 7 t) — q(j ± l,t) for the backward and forward shifts along the 
lattice, respectively. With these notations the equations of motion ^ become 

%r(j,t) = (S + -2 + S-)V'(r{j,t)). (11) 
One can write this equation as a first-order Hamiltonian system of the form 

d t u= ( x _ 5 _ 5+ _1 ) DH ( U )> ( 12 ) 

where u = (r,p), H(u) = £. £Z (fMi,*) 2 + V(r(j,t))), and DH(u) = (V'(r),p) 
denotes the functional gradient of H. By making the travelling wave ansatz r(j,t) = 
r c {j — ct) equation (12) becomes 

c 2 r''(x) = (S+-2 + S-)V\r c {x)). (13) 

Since we aim to consider spatially periodic solutions, the function spaces 
appropriate for our setting will be the periodic Sobolev spaces 

H$ L (R):={d>eV f 2L :U\\ 2 HL =2L £ (l + (^)*)|£( m )f < oo}, 

m— — oo 

where V' 2 i, is the set of 2L-periodic distributions, i.e. the set of all continuous linear 
functionals from V 2 l = C|l (^ ne se * °f a ^ smooth 2L-periodic functions from R into 
C) into C. The Fourier series of r c £ H\ L {M) and its Fourier coefficients r c are given 

by 

r c (x) = ^2 r c (m)e i!I ^ x , f c (m) = ^ / r c (a0e _i Tr x d:E.(14) 
The space if^^R) is a Hilbert space with the inner product 



m— — oo 

Remark 2.1. Formally, as first noted by Zabusky and Kruskal [50] the KdV equation 
arises from (13) as follows. One assumes that there exist solutions satisfying the 
multiscale ansatz 

r c (z)=e 2 <f(e(x-ct)), | = 1 + (15) 

which describes waves travelling at approximately the speed of sound while showing 
a negligible temporal and spatial change in their shape. Then by Taylor expanding 
differences and neglecting terms of order 0(e 8 ) in equation (13) one obtains a KdV 
equation - arising as the coefficients of the e e -terms - of the form 

-*" + 6 ^ (<I>2) " + $ "" = - (16) 

From comparison of the e^-terms one obtains the relation c s = ^V" (Q) for the speed 
of sound. 
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The ansatz ( 15 ) can be justified rigorously via a renormalization appraoch adapted 
from [T2] (see Section [3]) , spectral theory for the Schrodinger operator with KdV 
cnoidal wave potential (see Section [4] and Appendix A), and an implicit function 
theorem argument (see Section [5]) . This leads to the following rigorous persistence 
result for cnoidal waves when deforming the limiting KdV equation back into the 
lattice equation. 

Theorem 2.2. A ssum e that the nearest neighbour interaction potential V satisfies 
the assumptions (10), and let c s = v/V"(0) denote the sonic wave velocity. 

€ (0,1) and Lq G R + such that CKdv{ka,Lo) = 
1, where K(k) denotes the complete elliptic integral of 



Fix the constants ko 

±K 2 (koWl-k* + k*/L 2 

the first kind, i.e. K(k) = J^ 2 ds/\/l — k 2 sin s. Then the following statements hold: 

(a) (existence and local uniqueness) Given 6 > sufficiently small, there exists e > 
such that, fore € (0, £o) andc 2 /c 2 = (l + e 2 /12) ; the governing equation (13) for 
the periodic travelling wave profile admits a unique solution r c in the set 

{r G H^ Lo (R)\r even, H^r^- 1 .) - ^*°' io) || H x < 5}, 



where ^ l ko ' L °" > denotes the KdV cnoidal wave profile with speed CKdv (ko, Lq) = 1 
that solves the integrated KdV travelling wave equation 

V'" 1 



- $ + $' 



6 



Explicitly, $^ fc ° ,i °'' is given by 



V"(0) 
V"(0) K 2 (k 



2fcg 



k 2 



k 



V"'(0) 



T 2 



(17) 



+ Ken 2 




Here cn denotes one of the Jacobian elliptic functions (see Section^for its precise 
definition). 

(b) (asymptotic shape) The solution r c from (a) satisfies the estimate 



< Ce 1 



(19) 



H 



where C is independent of e. 
(c) (smoothness) The mapping c ^ r c from (c s ,oo) 



into FL\ Lq 



is C 1 



Estimate ( 19 ) means that the wave profile r c that solves ( 13 ) has a characteristic 
period of order \ je and amplitude of order e 2 . 



,(18) 



3. Renormalization and the continuum limit 



The formal argument in Remark |2.1| which relates the lattice equation to the KdV 
equation is not mathematically satisfactory. This is because it involves uncontrolled 
truncation of a Taylor expansion of a difference operator into (more and more 
unbounded) difference operators. To overcome this problem, Friesecke and Pego [12] 
introduced - in the context of spatially localized waves - a renormalization approach 



which recasts both the lattice travelling wave equation (13) and the limiting KdV 



travelling wave equation in a form involving only bounded operators (see ( 20 ) below) 



The ensuing lattice Fourier multiplication operator governing lattice waves and its 
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limiting continuum counterpart can then be shown, via a careful analysis of the 
location of their complex poles, to be rigorously close in an appropriate operator 
norm. As we show in this section, this approach can be adapted to the periodic 
setting. However, as we will see in Section|4j the limiting KdV equation ( 17) becomes 
more subtle in the periodic setting, admitting a two-parameter family of solutions in 
place of the one-parameter soliton family. 

We now derive the renormalized form of ( 13 ) and its small-amplitude long-wave 
limit, by adapting the analogous steps in [T^] to the periodic case. 

Let us make the assumptions ( 10 1 on the nearest-neighbour potential V. We 
isolate the nonlinear from the linear part of the restoring force by writing 

V'(r) = V"(0)r + N(r), N(r) = ^V'"(0)r 2 (l + j](r)), 

where rj(r) £ C 1 with 77(0) = 0. Let us consider functions r c in the periodic Sobolev 
space H\ L (M.) and their corresponding Fourier transform ( 14 ) . Then ( 13 1 turns into 
the equation 

( c2( "^) 2 _ 4 y"(0) sin 2 (~))r c (m) = 4sin 2 (^)iv^)(m). 

Since c 2 (=f f - 4F"(0) sin 2 (f^) > for c 2 > 1 and m £ 1 \ {0}, this equation can 
be transformed into the fixed point equation 

r c = P N{r c ), (20) 
where P denotes the pseudo-differential operator Pr(m) — p(m)f(m) with symbol 

, % 4sin 2 (^) sinc 2 (^) , 

n(m) = v 2L I = v 2L 1 (21) 

c 2(»hl) 2 _ AV»{Q) sin 2 C 2-F"(0)sinc 2 (ff)' 

Here sinc(a;) = sin(a;)/x denotes the usual sine function. 
Let us introduce a small parameter e > by setting 

c = VVW) (l + , with c Kdv = ^ Vl-k^ + k^, (22) 

where k 2 £ (0, 1) and L £ R + are fixed real constants. By rescaling variables via 

&*M(x) = \r c (-) (23) 



the renormalization equation (20) turns into an (e-dependent) fixed point equation for 

$(e,fc,£). 

$ = P (<r) iV( £ )($), (24) 

where N^(<^) = \V" '(0)$ 2 (l + rj(e 2 ^)) and the operator P^ has the symbol 

/ \ . e 2 sinc 2 (e^) , . 

y V ' c 2 -^"(0)sinc 2 (eff) 1 ' 

Next, we recall from [12] how the system ( [24] ) , (25) formally reduces to KdV as 
e — > 0. In the small e regime the wave speed scaling (22) implies 

-2, 



c 2 = V"(0)[ l + ^^l+O^ 4 ). (26) 



= V"(0) (l + ^)+0(e 4 ) 
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From this and the Taylor expansion sinc 2 (a;) = 1 — x 2 /3 + 0(x 4 ) one obtains for the 
pointwise limit of of the symbol as e —> 0: 

^w-^w.^—i-^. ,27) 

The nonlinearity satisfies lim e ^ A^ e ' (<&) = 5^"'(0)$ 2 . Hence as e -> the fixed 
point equation ( 24 ) converges to the equation 

$ = p(0) JV W($) = 6^±(c KdV - a 2 )- 1 * 2 . (28) 

By applying the operator (cxdv — d 2 ) to both sides, we see that this equation is 
equivalent to the integrated KdV equation 

V"'(Q) 

-Ci^ + e^-^^ + ^'-O, (29) 

which stems from integrating the usual KdV travelling wave equation 

V"'(0) 

-c KdV & + 12y^$& + &" = 0, (30) 

with vanishing integration constant. 

A key advantage of the above analysis over the formal asymptotics in Section [2] 
is the fact, shown in |12) . that the convergence in (27) is in fact uniform. This will 
allow us to show that the multipliers P^ and P^^ are close in operator norm on 
the periodic Sobolev space H\ L (J£) introduced in the previous section. It is useful to 
introduce the notation 

p (e) {m)=:p (e) { ^ p (0) (m)=:p -(0) ( !^ ); (31) 

the L-independent functions p^ and p^ coincide with those considered in [12] . 

Lemma 3.1. (Fourier multiplier estimate) Fix a point (fcoj^o) € (0,1) x R + such 
that CKdvi^Oj Lq) = 1. Then there exists a constant C > such that for all sufficiently 
small e > and c 2 = V"(0)(l + e 2 /12) the Fourier multipliers defined in (25), (21). 
pi] ) satisfy 

a) sup|p (£) (s)-p (0) (s)| <Ce 2 . 

b) \\P^~P m \\ £{HiLg) <Ce 2 , 

where C(X) denotes the Banach space of bounded linear operators from the Banach 
space X into itself. 

Proof, a) was proved in p~2j Lemma 3.1], by careful analysis of the location of the 
poles of p( £ ) in the complex plane. To show b), note that, for any operator on H^QSL) 
of form Ar(m) — a(rmr / L)f(m) , 



lML(H - = Jil 2L EmeZ (l + (W^) 2 )l^)l 2 = s!l /L WS)I 



Applying this with A = P^ - p(°), a = p^ - p(°) yields 



l^ (£) -^ (0) ll£ ( ^ ) = sup r 



Estimating the above supremum by that over s G R and applying a) yields the 
assertion. □ 



Cnoidal Waves on Fermi- Pasta-Ulam Lattices 



10 



4. Periodic KdV travelling wave solutions 

In this section we explicitly describe the periodic travelling wave solutions of the 



(integrated) KdV equation (29), see also OHO]- Here a main difference from the 



solitary wave case is present. Mathematically, there arises a non-zero integration 



constant after multiplying the KdV equation (17) by <&' and integrating. This forces 
us to deal with a 2-parameter system (for instance represented physically by the speed 
and period of the wavetrain) in contrast to the soliton case, which is fully described by 
one parameter corresponding to the velocity (or equivalently the amplitude or width) 
of the single soliton. The class of periodic solutions contains soliton solutions as a 
degenerate limit, namely infinite period. 



By multiplying the (integrated) KdV equation ( 29 ) with $' and integrating once 
more, we get 

^2_^"'(0) ( , 3 . V"(0) _ , , V"(0) 



(*T = 4^ (-^ + W^fK^ + j^S. ) ■ (32) 
Let us denote the polynomial in the bracket on the right hand side by F(&). For 



finding periodic solutions of (32) it is useful to consider the roots of F(<&). The 
function can by factorized by 

F($) = - £i)($ - £ 2 )($ - Ek), 

where E\, E 2 , and E 3 denoted the three roots. It turns out that real periodic solutions 
occur if the three zeros of the polynomial F($>) are real and distinct, i.e. E\ < E 2 < E 3 . 
The dependence of the zeros on the parameters ck<iv and B$> is thus given by 

V"(0) 

E l +E 2 + E 3 = W m^ c Kdv, (33) 

E 1 E 2 +E 2 E 3 + E 1 E 3 = Q 1 (34) 
V"(0) 

E\E 2 E 3 = ^5*. (35) 

The real and bounded solutions are contained in the region [E 2 ,E 3 ], whereas in the 
region (— 00, Ei] the solutions are unbounded. Introduce the normalized variable 



p = $/E 3 , then (32) becomes 



v"'((i) 

(pT = -^y^^P - pMp - (»)(p - !). ( 36 ) 

where we have set pi — Ei/E 3 for i = 2,3. The new variable p takes values in the 
interval [p 2 , 1]. Assume that a maximum value of p is at z = which can always be 
achieved by a coordinate translation. Since the critical points of p £ [p 2 , 1] are at the 
boundary of the interval we conclude p(0) = 1. Moreover, we have p"(p 2 ) > and 
p"(l) < 0. Hence a bounded solution p oscillates between the values p 2 and 1. 

In order to construct such a solution let us again make a change of variables by 
defining ip implicitly via 

p = 1 + (p 2 - 1) sin 2 V> 



where tp is continuous with -0(0) = 0. With this definition (36) becomes 



Cnoidal Waves on Fermi- Pasta-Ulam Lattices 



11 



We want to use the standard notation in this context and therefore define the 
parameters 



,2 x V"'(0) 



which fulfill k 2 € (0, 1) and A > 0. Using this notation (37) can be written in the form 

O') 2 = A(l - fc 2 sm 2 ^), 
which can be solved implicitly by integration 

ds 

mk 2 ):= =VAe + £o, (38) 

Jo yl-k 2 sin s 

where £o is a final arbitrary constant of integration. Here k 2 ) is called the 
standard elliptic integral of the first kind and k 2 the elliptic modulus. The Jacobian 
elliptic function sn is defined as the inverse of the function ip > /(?/>; fc 2 ) via 
sin?/> = sn(-\/A£ + £ ; k 2 ). Another basic elliptic function, the cnoidal function cn, 
is defined in terms of sn via cn(z; k 2 ) = yjl — sn 2 (z; k 2 ). Thus 

P = 1 + (p* - l)sn 2 (7Ae; k 2 ) = p 2 + (1 - P2 )cn 2 (\/A£ + £ ; k 2 ). 



Hence by recasting variables appropriately the solution of ([29]) is finally given by the 
so-called cnoidal wave 



S(0 = E 2 + (E 3 - E ^ 2 (\jy^( E ^ -E^ + ^k 2 ), (39) 

where the travelling wave coordinate is 

4V" (0) 

and the elliptic modulus is given by 

,2 E% — E-2 
k 



Ez — Ei 

The solution (39) is periodic if k 2 £ (0, 1) with period 

2L = 2,/V{G)K(k) 

^Vi'mEz-E^ 

where 

r /2 rfs 

K(k) := /(tt/2; k 2 ) = / ' 

Jo yl — k J sin s 

is called the complete elliptic integral of the first kind. 

It is worthwhile to consider two limiting cases of these KdV cnoidal wave solutions. 

(i) The "most nonlinear" limit: In the case E 2 — > E\ (corresponding to k 2 — > 1) we 
have K(k) — > +oo and cn(z;k 2 ) —> sechz. Thus in this case the solution of (16) 
with £q = corresponds to a single soliton solution 



$(0 = E 1 + (E 3 - E^sech 2 ^J^^(E 3 - E x ){x - c Kdv t) ) , 



'V"'(0) 
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where c Kdv = W"'(0)(2E 1 + E 3 )/V"(0). Due to the invariance of the KdV 
equation under Galilean transformations, i.e. transformations of the form 

/ V"'(0) 



V"(0) 

this solution is equivalent to the well-known standard form 

*« = H(f«Mf<*-w)) 2 , m 

where p = W"{0){E 3 - £?i)/V"(0). 

(ii) The "linear" limit: In the case E3 — > E2 which corresponds to the limit 
k 2 —> (amplitude tending to zero) we have K(k) —> tt/2 and cn(z;& 2 ) —> cosz. 



Therefore the limiting behaviour of ( 16 1 with £0 = as E 3 —> E2 is 



*(£) - ^+ <^> c„ 8 (f^A - E Ml - CKivt) ), (41) 

where c K dv = W"(0)(E! + 2E 2 )/V"(0). Here we made use of cos 2 (x) = 
(cos(2ir) + l)/2. Thus in leading orders of this small amplitude limit the cnoidal 
wave describes a linear wave. 

From the considerations above we deduce that the periodic solutions of the 



type ( 39 1 can be considered an intermediate form between linear waves and highly 



nonlinear solitons. 

By solving the KdV travelling wave equation (30) via integration one has in 
general to deal with a 3-parameter system, stemming from the cubic polynomial F(<f>). 
Thus the KdV solutions constructed in this way depend on 3 parameters, reflected 
by the roots {E\, E2, E3}, the integration constants and velocity {cxdv-, ^4*) Bq,}, 
or elliptic modulus, period and velocity {k 2 , 2L, CKdv}- Due to the choice of the 



renormalization ansatz &23h the governing equations ( 13 ) for the string are related to 



an integrated KdV equation (29) with zero integration constant, i.e. ^4$ = 0. Hence 
the 3-parameter reduces to a 2-parameter system. We aim to admit periodic KdV 
travelling waves to our considerations; for that purpose it is essential that the second 
integration constant is non-zero. Otherwise the system would include only one 
free parameter giving rise to soliton solutions. 



The cnoidal wave solution ( 39 1 is given in terms of the roots E\ , E%, and 
E 3 = E 3 (Ei, E 2 ). In the following we derive an equivalent expression in the parameters 
k, L, and cxdv = cxdv{k, L). Since the latter model directly the "physical" quantities 
like the shape, period, and velocity of the wavetrain this will help us to get a less 
abstract picture of the solutions we are dealing with. Besides it enables us to directly 



relate our spectral analysis results derived in Section 5.1 to studies on the Lame 
equation, cf. |Appendix A| 

Let us start by making the ansatz <&(£) = A+Bcn 2 (D£ > ; k 2 ) with £ = x—Ct. Insert 
expressic 
dn 2 x = ] 
of equations 



this expression into (29) and use the relations cn'x = — cnic sna; dna;, sn 2 x + cn 2 a; = 1, 
and dn 2 x — 1 — k 2 + k 2 cn 2 x. Then comparing coefficients yields the following system 



6 L F ^A 2 -CA + 2(l-k 2 )BD 2 =0, (42) 
V"'(0) 

c - l2 Y^ A + A{l ' 2k2)D2 = ^ (43) 
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0. 
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(44) 



Since we assume $ to be periodic with period 2L and cn 2 is a 2if(/c)-periodic function 
we have D = K(k)/L. Hence from (44) we get 

V"(0) K 2 (k)k 2 



B = 



V"'(0) i 2 



Moreover, from (43 1 one obtains 
V"'(0) 



C = 12 



V"(0) 



A _4(l-2fc 2 )^P. 



Finally, inserting these formulas into ( 44 ) gives the following quadratic equation 

2 



A 2 - 2(1 - 2k 2 



V"(0) K 2 (k) 



V"'(0) 3L 2 
which has the positive solution 

V"(0) K 2 {k) 



A-(k 2 - fc 4 ) 



V"(0) V K ( k ) 



V"'(0) 



3L 4 



.4 



(l - 2k 2 + v 7 ! - k 2 



k 4 



V"'(0) 3L 2 

Thus we can write the cnoidal wave solution ( 39 ) in the convenient form 
V"(0) K 2 (k) ( 1 - 2fc 2 + Vl - k 2 + k 4 



V"'(0) L 2 



+ k 2 crt 



( K(k) 



£;fc 2 ,(45) 



where £ 



c Kdvt is the travelling wave coordinate with the velocity 



CKdV 



AK 2 (k) 
L 2 



Vl-k 2 + k 4 . 



Note that for k £ (0, 1) we have cxdv € (tt /L , oo). The representation (45) of the 
cnoidal wave solution is indeed equivalent to ( 39 ) if we choose 



E* 



V"{0)K 2 (k) /r -^ r 
~ V"'(0) 3L 2 [1 2k +V1 k 
V"(0) K 2 {k) 



fc 4 ), 



(l + fc 2 + Vl-^ + fc 4 )- 



y'"(0) 3L 2 

In Figure [2] cnoidal wave solutions of the form ( [45] ) for different parameters A 2 
and c^dv are shown, that is they are normalized to have the same period. 

Remark 4.1. Note that an equivalent re nor malization approach would be the 
following: replace the renormalization ansatz (23) by (f>^ e ' k ' L \x) + E% = e~ 2 r c (e~ 1 x), 
where E 2 7^ denotes an arbitrary constant. In the limit e — > this leads to an 
integrated KdV equation for <f> given by 



-CKdv4> + 6 



V"'(0) 



A, 



0. 



V"(0) r 

where A^ 7^ 0. Indeed, just set $ = + E 2 in When integrating this equation 

once more assume the second integration constant B^ to vanish. Then one of the zeros 
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Figure 2. Cnoidal wave solutions | |45| l of | |29| l normalized to have the same 
period. 2L = 2tt. We have set V"(0) = V"'(0) = 1. The three curves represent 
the solutions corresponding to the different values for the elliptic modulus k 2 = 0.6 
(dotted), k 2 = 0.9 (dashed), and k 2 = 0.999 (solid). 



ei < e 2 < e3 of the corresponding cubic polynomial has to vanish, say e 2 = 0. The 
cnoidal wave solution is given by 

<H£) - e 3 cn 2 [y/V"(0)(e 3 - ei )/V"'(0)e, fc 2 ) , with £ = x - c Kdv t, 

where the velocity is ck<iv = 4V A '"(0)(ei + e%) /V" {0) . If we now add the constant 
E%, that is if we make the transformation <p — > <p + E2, one has to shift the velocity 
simultaneously, C-KdV cxdv + 12V r "'(0)i?2/V r "(0) ; to obtain an equivalent solution 
( Galilean invariance) . Indeed, this yields the solution $ given by (39) with Ei = e.;+i?2 
fori = 1,2,3. 



5. Existence and shape near the continuum limit 



In this section we prove our main result, Theorem |2.2| The technical ingredients of 
the proof are 

(i) the renormalized form of the lattice travelling wave equation and associated lattice 
Fourier multiplier estimate given in Section [3j 



(ii) spectral information about the linearization of (20 1 at the KdV cnoidal wave 



(results from direct scattering theory for Schrodinger operators with KdV cnoidal 
wave-potential), and 

(iii) a quantified version of the standard implicit function theorem borrowed from |12) . 

We start by fixing the wave speed ck<w of the KdV travelling wave profile by fixing 
the two parameters k and L of our system. We do so by choosing a point k Q £ (0, 1) 

1. That means, we 



4if 2 (fc )\/l - ^0 + K sucn tnat c Kdv{ko, Lq) 



and setting Lq 
fix the parameters k and L in (|23| and (|26|) 
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In Figure [T] cnoidal wave solutions of the form (45) for different parameters k 2 are 
shown. By setting L? = AK 2 (k)\/l — k 2 + k A they are normalized to have the velocity 

C-KdV = 1- 

Next, we introduce 

= p(°'iV (0) ($), F($) = p( £ 'iV (£) ($) - P (0) A^ (0) ($). 

The map F is smooth on FL\ La and it has a fixed point <&V so,Lo ' given by (45). The 
derivative of F at $( fc °' Lo) is the operator L = DF(<f>{ k °' L " ] ) defined by 

V"'(0) $^°' Lo) 

We will solve the equation $ — F($) — i^(<I>) = by the inverse function theorem, using 



that (i) F is small in a suitable sense for small e (using the estimate of Lemma 3.1 ) 



and that (ii) the linearization of the map <I> - - at the fixed point ${ k °' 

for e = is I — L, which has a bounded invers e on the space of even 2Lo _ periodic 
functions E\ Lq = {$ £ H\ Lo \<& even}, cf. Lemma 



5.1 



5.1. Spectral analysis 

To apply the implicit function theorem we have to establish spectral properties of the 
operator = pW N(°\&) linearized at a KdV cnoidal wave solution ^ k ' L \ For 

our purpose it suffices to show that the eigenvalue A = 1 is simple with corresponding 
odd 2L-periodic eigenfunction ^$[ k ' L \ such that / — DF ($[ k ' L ^) has a bounded 
inverse on the space of even 2L-periodic functions. 

It turns out that the spectral problem for L = DF(&[ k ' LS> ) is connected with the 
scattering theory for Schrodinger operators with cnoidal wave potential. There exists 
a large literature on this eigenvalue equation, known as the Lame equation. Therefore, 
an alternative way to obtain the below result via the theory of the Lame equation is 
given in | Appendix A 



Lemma 5.1. Fix a point (ko,Lo) £ (0,1) x R + such that CKdv(ko, Lq) = 1. Then 



the operator L = D_F($^ fc °' Lo '') given by (4.6) is a compact self-adjoint operator 
on FL\ L (K). The eigenvalues X% = 1 and A2 = 2 are simple with corresponding 

eigenf unctions — (£) and V*2(£) = $^ fc °' L °' (£), respectively. Hence 

the operator I — L has a bounded inverse on the space of even 2Lo~periodic functions 
E \l q = {<l>£HL\<t> even}. 



Proof. That L is a self- adjoint operator follows from 

+00 



(Li>,4>) H i =2L (1 + (m7r/i ) 2 )^(m)?(m) 

771= — OO 

=2io3 ) g K°' Lo v)M^= i ^r o (^ feo - io) ^)(,)^ 

^ ' m— — 00 ' ' ' l; 



The compactness of L is a consequence of the fact that L is a bounded operator 
from FL\ La into H 2 Lg and the Rellich compactness theorem, which guarantees that 
H\ L is compactly embedded into H\ La ■ 

Concerning the spectral properties of L we observe the following: since L is a 
compact self-adjoint operator all eigenvalues are real and semi-simple. Indeed, all 
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eigenvalues are simple: note that the spectral problem Lip — \ip where L is given by 
( 46 1 is equivalent to the second order differential equation 

12 V"'(0) 



CKdVlp 



A V"(0) 

which can n ot have two independent bounded solutions. Moreover, if A 



(47) 



1 then a 



d ^(feo.io) 



since $ 



(ka,L ) 



solution of (47) is if) = J|$ 

(301. If A = 2 then a solution of (|47j is ip = $i' 
integrated KdV equation ( 29 1 . 



(ko,L, 



is a solution of the KdV equation 
) ' ) since $^ ' L °' is a solution of the 



□ 



5.2. A quantified version of the standard implicit function theorem 

The following result will be borrowed from [T^], where it has been applied to the 
solitary wave case. 

Lemma 5.2 (Quantitative version of the implicit function theorem, Lemma A.l). 
Suppose F and F are C 1 maps from a ball in a Banach space E into E . Suppose F 
has a fixed point $i 7 i.e. F(<f>±) — $i. Let L = DF (<f>i) denote the linearization of F 
at $x- Then assume that the operator L — L is invertible with \(I — L)~ 1 \ < Co < oo. 
Assume that S > is sufficiently small such that 

\DF($) - DF($i)| < Cx < Cq 1 for \<f> - <f>x\ < 5. 

Suppose that F is small in the sense that for some positive constants Ci and 
satisfying Cq(C\ + C'2) < 9 < 1, we have 

|F($i)| < 6(1 -8)/C and \DF(<f>)\<C 2 for |$-$i|<<5. 

Then the equation $ = F(<$>) + F(<f>) has a unique fixed point satisfying |$ — $i| < 5, 
and moreover 

|$-*i|<(l-fl) -1 Cb|F(*i)|<*. (48) 



5. 3. Proof of the main result 



Here we finally give the proof of our main result Theorem 2.2 by combining the Fourier 
multiplier estimate from Lemma |3.1[ the spectral information gathered in Lemma |5.1[ 
and the above version of the implicit function theorem. 



Proof of Theorem \2.2\ We first prove statements (a) and (b) . Introduce the operators 
F($) = p(°)Ar(°)($), F($) = P (e) iV^($) - PWJV(°)($), 



where the symbols of the operators P^ and P^ are given by (25) and (27), 
respectively. Moreover, recall that 



We now want to solve the fixed point equation |24|) , which is equivalent to the equation 
$ = F(&) + P(<I>), by the implicit function theorem. Hence we prove that on the space 
E\ L (K) of even functions in FL\ L (R) the operators F and F satisfy the criteria of 
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Lemma [572] That F fulfills the required properties has been shown in Lemma [5. 1| It 
remains to establish the estimates on F and DF. For that purpose write 

= p( £ )jv( £ >($) - p(°)ivW($) 

= p( £ )(iv( £ )($) - at(°)($)) + (p( £ ) - p(°))iv( )($) 

= ^)^1$^( £ 2$) + ( P W _ p(0))Z^)$2 

Now choose 5 > sufficiently small such that the hypotheses of the Lemma |5.2| are 
guaranteed. Since rj € C 1 with 7/(0) = we have ||?7(e 2< I>)|| < Ce 2 . Moreover, using 
Lemma 3.1 b) we deduce 

\\F\\c(H ko ) < Ce\ ll^l| £ (^ o ) < Ce\ 



where C is a constant indepe nde nt of e. Then Lemma |5 . 2 1 implies that there exists a 
unique solution $( e ' fc °' L °) f (24) satisfying 

|| $ (e,fc ,io) _ $( fe o,io)|| < 



where $( fco,L °) is the fixed point solution of $ = P($) and given by (45 ) with ko and £*o 
hxed such that CKdv{k ,L ) = 1. With r c := s 2 ^ e ' ko,L °^(£.) this proves statements 
(a) and (b) of the theorem. 

For proving the statement (c), we fix e but reintroduce the parameter CKdv = 
AK 2 (k)^J\ — k 2 + k 4 /L 2 by assuming the parameters k £ (0, 1) and L € M+ not to be 
fixed. That is, we take c 2 = V^"(0)(1 + e 2 ck<iv /^)- Then we define the map 

G($,fc,L) = P (e) iV (£) ($), G : E\ L x (0,1) x M+ — > E\ L . 

Now observe that {k, L) i-> cxdv and thus (fc, i) i— > c is C . Thus the mapping 

(fc, L) ^ p = sinc 2 ( • ^)/(c 2 - ^"(0)sinc 2 ( • ^)) 

is smooth from (0, 1) x K + into L°°(]]i). From this it follows that the map G is C 1 . 
For (kQ,Lo) fixed such that cxdv = 1 we have shown that there is a fixed point 
^(e.feo^o) anc j f rom existence proof it follows that I — D$G($> < f' ko ' L °\ ko, Lq) is an 
isomorphism on E\ Lo . 

From the standard implicit function theorem [51] Theorem 4B] we can thus deduce 
the following: for some small enough neighbourhood of (k a ,L ) g (0,1) x M + there 
exists a C^-curve (k, L) h-> $( £ - fc < L ) of solutions of Q such that $( e <^ L o) _ ^,k ,L )^ 
Thus the mapping (fc, i) h-> c^ v $( e,fc ' L )(c^y .) is C 1 using again that (fc, i) H> c^d^ 
is C 1 . Note that c^ y $(^ fe ^) (c" 1 /^.) is a solution of $ = p( £ i)A^i)($) w jth £l = 
^JCKdv£ and c 2 = V" (0)(1 + e 2 /12). Moreover, if is sufficiently small and (fc, L) in 
a small enough neighbourhood of (k ,L ) then \\c K l dv <S> {e > k > L) (c^Jy.) - <f>[ £ ' k °> L ^ || < 5 

and so c^& e ' k,h) ( c ~Kdv-) = <S> (El ' k ' L) - It follows that (k,L) i— > and 
hence c = (V"(0)(1 + e 2 /12)) 1 / 2 ^ r c = ef $( El ' fe ' i )(e 1 .) are C 1 maps into P^. □ 
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Appendix A. Results from the theory on the Lame equation 



The spectral problem for L defined by ( 46 1 which reads 



h x = 



L^(0=AV(0 on Hi L (R), (A.l) 
is equivalent to the following problem 

H^(y) = h^(y) on H^ K (R), (A.2) 
where the transformed coordinate, operator, and eigenvalue are 

y = *t, (A.3) 

H^= ~d 2 -^fc 2 cn 2 (y;fc 2 ), (A.4) 

i(4 - 8fc 2 ) + 4^ - l) v 7 ! - fc 2 + A: 4 . (A.5) 
Using the relation cn 2 + sn = 1 this spectral problem transforms into 

H^(y) = h^(y) on H$ K (R), (A.6) 
where the operator and eigenvalue are 

-a 2 + ^fc 2 sn 2 (y;fc 2 ), (A.7) 

hL = h x cn + ^k 2 . (A.8) 

The general eigenvalue problem for the operator H^ n is connected to the Jacobian 
form of the Lame equation, namely 

- ip"(y) + n(n + l)fc 2 sn 2 (y; k 2 )^j{y) = Eip(y). (A.9) 

For n € N the solutions of this equation are meromorphic in the whole complex plane, 
while for n other than an integer they are multi-valued. We will use results from the 
theory on the Lame equation [MJ [5] combined with well-known facts from Floquet 
theory [51]. 

We start by recalling the basic facts from Floquet theory. Let us consider the 



spectral problem (A.9) together with the boundary condition 

4>(y + 2K) =e i2KK i;(y) := a^(y), 

where k is referred to as the momentum and <j the Floquet multiplier. For any 
eel, i.e. \a\ = 1, this defines a self-adjoint boundary value problem. The spectrum 
of the 'periodic (a — 1) and semi-periodic (a = —1) problem consists of infinitely 
many discrete eigenvalues which only accumulate at +oo and the corresponding 
eigenfunctions are 2K- and 4if -periodic, respectively. Let Eq < E^ < E~£ < ■ • ■ 
denote the eigenvalues and i/Jq , ip^ , ip^ , ■ ■ ■ the 2 if -periodic eigenfunctions of the 
periodic problem. Similarly, Eq < E^ < E% < ... and ipQ , ip± , ip% > • • ■ are the 
eigenvalues and 4i4T-periodic eigenfunctions of the semi-periodic problem. We have 
that 

(i) ipQ has no zeros and i/jq exactly one zero in [0, 2K(k)], 

(ii) tp2i-i anc ^ ^ti nave exactly 2i zeros in the interval [0, 2K), for i > 1, and 

(iii) 1P21-1 ano - V"2i nave exactly 2i — 1 zeros in the interval [0, 2K), for i > 1. 
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Then by oscillation theory we have 

Eq < Eq < E± <E+ <E+ <Ez <E^ < ... 

The bands 

(E+,Eq), (E^,E+), (E+,Ez),.., 

are the stable bands. For each value in these bands we have k€R and hence \o~\ = 1, 
i.e. only bounded but non-periodic solutions exist. Inbetween these bands we have 
gaps 

i E 0, E o)' ( E l> E l)> i E 2l E 2~)>--' 

which are called bands of instability. For values of E lying in these gaps only 
unbounded non-periodic solutions exist, since the multiplier \a\ ^ 1 and the 
momentum k can not be real. 

The solutions of the Lame equation (A. 9) are given by 

<i>± (») = n H %tv) j) exp ^ yz ^ ' 

for certain constants a±, . . . ,a n , see [49]. Here H denotes the Jacobian eta, 8 
the Jacobian theta, and Z the Jacobian zeta function [5]. At a band edge the 
eigenfunctions correspond to polynomials in the Jacobian elliptic functions cn, sn, 
dn and are called Lame polynomials. 

It is known [53], [Ml Theorem 7.8] that if and only if n is an integer, 2K- or 4K- 
periodic solutions can coexist, that is the spectrum has finite-band structure (since 
two band edges collide) with n + 1 bands (if n is non-negative) . The 2n + l band edges 
correspond to simple eigenvalues of the periodic or semi-periodic problem. 

Let us consider the periodic problem 

H^{y) = H{y), V(0) = V(2#), V'(0) = i//(2K), (A.10) 

where H™ n equals the operator given by ( |A.7[ ) if we set 12/A = n(n + 1). Ince [23] 
studied the Lame equation for n G R, hence not necessarily an integer. From his 
results it follows that there exist 2K -periodic eigenfunctions of Hf n with corresponding 
eigenvalue h of the specific form (A.8), that is in terms of n 



^ (n(n + 1)(1 + k 2 ) + 4(n(n + 1) - 3) v / l-fc 2 + fc 4 ) , (A.ll) 



h" = 
sn 3 

if and only if n = 2 (A = 2) or n = 3 (A = 1). Only in these cases Ince's continued 
fraction representations for h terminate and give rise to solutions given by the finite 
form (A.ll |. In cither cases, n = 2 and n — 3, the concerned eigenvalues h 2 sn and 
hg n are simple and correspond to a band edge of the finite band spectrum of H 2 n and 
Hg n , respectively. Note that an eigenvalue A of the spectral problem (A.l) is simple 



if and only if the corresponding eigenvalue h^ n of (A. 6) is simple. 
Case A = 1: 



In this case the spectrum of the operator H^ n consists of the four spectral bands 

(E+ ,£?„), (£V , E+ ) , (£ 2 - , E+ ) , (£ 3 - , oo ) , 

where the E^'s are the eigenvalues corresponding to the periodic problem, and the 
i?~'s the ones correspondin g to t he semi-periodic problem. Hence the first three 
eigenvalues E^E^E^ of (A.10) are simple and all other eigenvalues E^ with 
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Case n = 3 (A = 1) 


eigenvalue 


eigenfunction 


period 


parity 


E+ = 2 + 


5fc 2 


- 2s/\ - 


k 2 -+ 


-4fc 4 


i>+(y) = dny[l - (E+ 


— k 2 )sn 2 y] 


2if 


even 


Eo=5 + 


2fc 2 


- 2V4- 


fc 2 + 


-fc 4 


</>o (y) = cnyl 1 - (^o 


— l)sn 2 y] 


4if 


even 


Ei = 5 + 


5fc 2 


- 2V4- 


Ik 2 


+ 4fc 4 


V^r (y) = 3sn y - (^r - 


- 1 - fc 2 )sn 3 y 


4if 


odd 


£/" = 4 + 


4fc 2 








^i" (y) = sn y cn y dn y 




2K 


odd 


= 2 + 


5fc 2 


+ 2Vl - 


fc 2 + 


-4fc 4 


i>+(y)=dny[l-(E+ 


— fc 2 )sn 2 y] 


2K 


even 


#2~ = 5 + 


2fc 2 


+ 2V4- 


fc 2 4 


-fc 4 


^2 (y) = cnyl 1 - ( E 2 


- l)sn 2 y] 


AK 


even 


E 3 " = 5 + 


5fc 2 


+ 2V4- 


7fc 2 


+ 4fc 4 


i>3 (y) = 3sn y - (-^J - 


- 1 — fc 2 )sn 3 y 


AK 


odd 



Table Al. The spectral quantities of the spectral problem of the operator 
Hf n = -&l + 12fc 2 sn 2 (y; k 2 ) with periodic or semi-periodic boundary conditions. 
The eigenvalues correspond to the band edges of the finite band spectrum. The 
associated eigenfunctions have period 2K or AK and are polynomials in the 
Jacobian elliptic functions sn(j/;£: 2 ), cn{y;k' 2 ), and dn(y;fc 2 ). 



i > 3 are double. From [24] or [2] (the eigenvalues and eigenfunctions are calculated 
explicitly for the cases n = 1, 2, 3) we get results summariz ed in T able Al 



The eigenvalue E^ = A + Ak 2 is indeed of the form (A.ll) for n — 3 and the 



corresponding eigenfunction has period 2K and is given by ipi(y) = sny cny dny. 
Thus tracing back definitions we get the following result on the spectral problem 
(A.l): the operator L has the simple eigenvalue A = 1 with the 2L-periodic odd 



eigenfunction ip(£) — sn(A£/L; k 2 )cn(K£/L; fc 2 )dn(i-T£/L; k 2 ). Indeed, this gives the 

To see this note that 



(cny)' 



-sny dny and 



same result as derived in Lemma 5.1 
hence ±§f ' L) oc ip(£). 

The eigenvalues of the periodic and semi-periodic spectral problem listed in 
depend on the elliptic modulus k 2 . This dependence is visualized in 

1 we have cny ss sechy. The finite band 



Table 
Figure 



Al 



Al 



In the limiting case k 2 
structure arising from the periodic potential turns into a purely discrete spectrum. 
The corresponding eigenfunctions are given by the associated Legendre polynomials 
in tanhy. In the case k 2 — > we have cny ss cosy and the Lame equation turns into 
Mathieu's equation. The eigenfunctions are then polynomials in the Mathieu cosine 
and sine function. 
Case A = 2: 

In this case the spectrum of the operator H 2 n consists of the three spectral bands 

(E+,E^),(E^,E+) 1 (E^,^), 

where the E^'s are the eigenvalues corresponding to the periodic problem, and 
the E~ 's the one correspo nding to the semi-periodic problem. Hence the first two 
eigenvalues Eq and E± of (A. 10 1 are simple and all other eigenvalues E^ with i > 2 
are double. From [23] or [5] we get the results summarized in Table A2 

The eigenvalue E^ = 2 + 2k 2 - 2^J\ - k 2 + fc 4 is indeed of the form < |A.ll ) 
for n — 2 and the corresponding eigenfunction has period 2K and is given by 
^o~(y) = 2 — £g~sn 2 y. Thus tracing back definitions we get the following result on 
the spectral problem (A.l): the operator L has the simple eigenvalue A = 2 with 



the 2L-periodic even eigenfunction ip(t;) = 2 — [1 + k 2 — \J\ — k 2 + fc 4 ]sn 2 (AT£/L; fc 2 ). 
Using the relation sn 2 y = 1 — cn 2 we get 3>i (£) oc verifying our result derived 

in Lemma 15. II 
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Figure Al. The finite band structure of the periodic operator H%„ . The band 
edges depend on the elliptic modulus k 2 £ (0, 1) as given in Table [Al] 



Case n = 2 (A = 2) 


eigenvalue 


eigenfunction 


period 


parity 


E+ = 2 + 2k 2 


- 2y/l - k 2 - 


h/c 4 


^+(y) = 2-E+sn 2 y 




2K 


even 


Eq = 1 + k 2 






(y) = cn y dn 2/ 




4K 


even 


Ei = 1 + 4/c 2 
Et = 4 + fc 2 






ipl (y) = sny dny 
ipi{y) = sny cny 




AK 
2K 


odd 
odd 


Ef = 2 + 2k 2 


+ 2V1 - k 2 - 


Vk A 


^+(y) = dny[l-(E+ - 


- k 2 )sn 2 y] 


2K 


even 



Table A2. The spectral quantities of the spectral problem of the operator 
H 2 n = —dy + 6fc 2 sn 2 (y; k 2 ) with periodic or semi-periodic boundary conditions. 
The eigenvalues correspond to the band edges of the finite band spectrum. The 
associated eigcnfunctions have period 2K or 4K and are polynomials in the 
Jacobian elliptic functions sn(y;k 2 ), cn(y;k 2 ), and dn(y;fc 2 ). 



Appendix B. Quasi-periodic KdV solutions 

In this appendix we give an overview of the beautiful theory of quasi-periodic KdV 
solutions. 

The inverse scattering transform developed by Gardner et al. |17j allows to solve 
the Cauchy problem for the KdV equation 

u t + 6uu x + u xxx — 0. (B.l) 

Unfortunately, this method is restricted to rapidly decreasing initial data. Due to the 
lack of methods for solving the inverse spectral problem for H with periodic potential 
it had not been immediately clear how to adopt the inverse spectral method to the 
periodic case. 

A breakthrough in this problem was achieved independently by Dubrovin [3 E] , 
Its and Matveev [26} [2Z| , and McKean and van Moerbeke [35] , generalizing an earlier 
result by Akhiezer |T] . They proposed a scheme to construct so-called finite-gap KdV 
solutions from given spectral data. 

To describe this method in a little more detail, consider the Lax operator with a 
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periodic potential, i.e., 

d 2 

H = — —r+u(x,0), u(x + 2L,0) = u(x,0). 
dx z 

Assume the absolutely continuous part of the spectrum of H has finite-band structure, 
i.e. it consists of N + 1 intervals [E^j-i, E%j] € R, j = 1, . . . , N + 1 with E2N+2 = 00. 
Then consider the Dirichlet problem 

Ei\) = p4, V(0,0)=V(2£,0) = 0. 

We refer to the spectrum of this problem as auxiliary spectrum consisting of a set of 
discrete points {fj,j(0,0)}°^ 1 . All but N points of this set lie in the infinite interval 
[E2N+1, 00) of the continuous spectrum. Each of the N remaining points is contained 
in one of the N spectral gaps. Finally, the solution u{x, t) is given by the trace formula 

2N+1 N 

u(x,t) = ^ E i -2^2n j (x,t). (B.2) 
i=i 3=1 

Here the functions fij(x,t), j = 1, ...,N denote the space-time evolution of the 
auxiliary spectrum. These quantities live on the Riemann surface associated with 
the curve 

2JV+1 

R : y 2 = - [] (B.3) 

j=i 

where each of them lies in a spectral gap, i.e. (ij(x,t) € {E2j, ^j+i), j = 1,...,N 
and oscillates between these band edges as x varies. More precisely, the dynamics in 
x and t is given by the Dubrovin equations 

where the prime and dot denote the derivative with respect to the space and time 
variable, resepctively. Solving these equations is a matter of applying the Jacobi 
inversion problem to the so-called Abel map, that is a map from the Riemann surface 
R to an iV-torus called the Jacobi variety. In the case N = 1, finding a periodic KdV 
solution is linked to the classical inversion problem for the elliptic integral, cf. Section[3] 



The expression (B.2 1 for the solutions of the KdV equation (B.l I, called finite-gap 



KdV solutions, is finally given by 

u(a; > t)=c+2jLlne(0, (B.5) 

G(£;B)= exp{27ri<k,0 + 7ri(k,Sk)} ) 
kez™ 

where (., .) denotes the scalar product in R N and c a constant. Moreover, B is 



the symmetric period matrix of the Riemann surface R of the curve (B.3) and 
£ = ,£jv) denotes the vector of the generalized phases £j = KjX — u>jt + 4>j, 

j = 1, . . . , N. The Kj's and ujj's are completely determined by the continuous and 
auxiliary spectrum and can be computed in terms of Abelian integrals. Note that the 
Riemann 0-function fulfills the property 

8(£ + m + Bn; B) = cxp {-27ri(n, £) - 7ri(n, Bn)} 9(£; B), 

for all n, m € Z N and symmetric N x N matrices B. 
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In the one-gap case (N — 1), the Riemann O-function reduces to one of the Jacobi 
^-functions which are linked to the Jacobi elliptic functions cn, sn, and dn [5]. The 



solution (B.5I is then given by 



u{x, t) = E 2 + (E 3 ~E 2 )cn 2 (^J Ez - El (x-2(E 1 +E 2 +E 3 )t) ; fc 2 ) , (B.6) 
where E\ < E 2 < E 3 and k 2 = (E 3 — E 2 )/(E 3 — E\) is the elliptic modulus. 



Equation (B.6 1 is the so-called cnoidal wave which describes a periodic travelling 



wave solution of (B.l I and has been already known to Korteweg and de Vries in 1895 
[50] . The infinite period limit, E\ — > E 2 (k 2 —*•!), gives the single soliton 




u(x,t)=PsectfU^(x-f3t)), (B.7) 
where we have set j3 — E 3 — E 2 and imposed the condition u — > as x — > ±oo. 



Note that any finite-gap solution of the form (B.5 1 is an Abelian functioijJJJ which 
constitutes a generalization of elliptic functions to TV variables. Thus, in the special 
case N = 1 the solution is elliptic and referred to as cnoidal wave or Lame potential. 

is periodic in x (since the initial condition is), but in general 



The solution (B.5 



is only quasi-period!^ in t. In general, however, any finite-gap KdV solution is quasi- 
periodic in both x and t, since their construction via algebro-geometric methods [3"Ill8| 
allows to choose the quantities {Ej} 2 ^ 1 £ K and {//j(0, 0)}^^ arbitrarily such that 
the Kj's and w^'s are in general incommensurate quantities. More precisely, these 
methods use results on Abelian functions and Riemann surfaces. A key role is played 
by the eigenfunctions of the spectral problem for the Lax operator H. They turn out 
to correspond to the two distinct branches of one single- valued function, the so-called 



Baker- Akhiezer function, defined on a two-sheeted Riemann surface given by (B.3l. 

Note that in the case of the the Toda lattice, where the interaction potential 
is given by V(r) = a(e~P r + fir — 1) and the system is integrable, there have been 
constructed exact cnoidal wave solutions |45j . For the construction of all quasi-periodic 
finite-gap solutions for the Toda lattice we refer to [33| or [TO] , 

Finally, we would like to point out two major points which emphasize the 
significance of finite-gap solutions for understanding KdV dynamics: 



(i) The finite-gap solution (B.5) degenerates to a multi-soliton solution if the band 
edges of the continuous spectrum and thus the branch points of the corresponding 
Riemann surface R collide, i.e. E 2 j — > E 2 j + i, j = 0, 1, . . . , N — 1. In this case 
the Riemann theta function reduces to a determinant with exponential entries. 
In this way the finite-gap KdV solutions ( |B.5[ ) corresponding to N gaps in the 
continuous spectrum constitute a natural (quasi-)periodic extension of the KdV 
iV-soliton solutions. 

(ii) The KdV equation has almost-periodi^^ infinite-gap solutions corresponding to 
an infinite number of gaps in the spectrum [32 . The set of all quasi-periodic 

|| An Abelian function is a meromorphic function in N complex variables with 2N independent 
periods, the so-called fundamental periods. It is usually defined as a function on an Abelian variety. 
1 A finite-gap solution is called periodic in x with period 2L if the wave numbers K i , . . . , ii"jv are 
commensurable, i.e. there exist TV integers m, . . . , rijv such that ILKj = n,-, j = 1, . . . , TV. Otherwise 
it is called quasi-periodic. 

+ A function / is called almost periodic in the sense of Bohr if for each e > there exists a 
relatively dense set of e-almost-periods, i.e. translations r with the property sup x (f(x + r), f(x)) < e. 
Alternatively, there is the Bochner criterion for almost periodicity: a continuous function / is almost 
periodic if and only if the family of functions H = {/ } = {f(x + h)}, — oo < h < oo, is compact in 
the set of continuous bounded functions. 
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finite-gap solutions (which also contains all periodic finite-gap solutions) is dense 
in the former set. This means, in particular, that any periodic KdV solution can 
be approximated by finite-gap KdV solutions These considerations concern 
the periodicity in the spatial variable, but what about the behaviour in time? 
It is known [SHI EI] that all KdV solutions which are periodic in x are almost 
periodic in t. 
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